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Abstract 

Based on the Keldysh nonequilibrium Green function (NGF) technique, a general formula for 
the current and transmission coefficient through a one dimensional lattice is derived without the 
consideration of electron-electron interactions. We obtain an analytical condition for perfect reso- 
nant transmission when the levels of sites are aligned, which depends on the parity of the number 
of sites. Localization-delocalization transition in a generic one dimensional disordered lattice is 
also analyzed, depending on the correlation among the hopping parameters and the strength of 
the coupling to reservoirs. The dependence of the number and lineshape of resonant transmission 
and linear conductance peaks on the structure parameters of the lattice is also given in several site 
cases. 
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I. INTRODUCTION 



Transport through a quantum dot lattice depends on the matching of electron levels in 
the various quantum dotsQ]. If the electron levels of individual dots are aligned, resonant 
tunneling occurs [2]. Electronic transport through these structures is usually investigated 
theoretically on the basis of a "classic" rate equation in the weak coupling limit P, 0], and 
a quantum rate equation in the strong and weak coupling limits |l L 2|. The splitting due 
to the interdot coupling has also been calculated by Matveev et al. 5| and Golden et al. f|. 
Some groups have recently investigated in detail the spectral and current properties of a 
coupled quantum dot pair following the Keldysh nonequilibrium Green's function (NGF) 
formalism ?]. Calculations on Phonon-assistedS] and photon- assisted transport in two 
coupled quantum dots have been recently conducted, in which an electron pump effect is 
found. Kondo physics and other correlation effects in the same system have also been 
considered In quantum dot systems, electron-electron interaction becomes important 



and Coulomb oscillation of linear conductance and even Kondo effect arise. 



.LllljJ Exper- 

n 

imentally, as early as 1990, Kouwenhoven et al. [13( performed a transport experiment in a 



sample consisting of fifteen dots. Some experiments have also been conducted to explore the 
ground-state properties of a double quantum dot[l4|. Very recently, Oosterkamp et al. [l^ 
investigated experimentally microwave spectroscopy of a two quantum dot molecule. 

Recentl y e lectronic transport in metallic chains of single atoms has attracted both 
theoreticaljlfil. \v\ and experimental interest, since atomic wires represent the ultimate 
limit of the miniaturization of electrical conductors. Although atomic wires are somewhat 
simple toy models in textbooks and literatures, they may be the best laboratory to test 
physical properties of one dimensional (ID) systems. The microscopic mechanism for the 
conductance quantization of atomic wires still keeps vague in spite of some endeavours !l 9]. 
First-principle calculations[16] and analytical deduction |Tq] demonstrated interesting os- 
cillation of conduction of the monoatomic wires with the number of atoms, which was 
observed experimentally in the conductance of mechanically controllable break junctions 

I L 

atomic contacts 18]. The parity effect in the conductance was also found in weakly dis- 
ordered quasi-lD tight-binding modelsj^ and dirty superconducting wires I21I1 . However, 
whether the interesting parity phenomena arising from single particle behavior |l7( or many- 
body correlation effects [221] need more clarifications. 
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To have a better understanding the transport properties of generic one diemnsional sys- 
tems and clarify conductance quantization in such a system, we consider electron transmis- 
sion through a one dimensional lattice where the sites allow for only the nearest-neighboring 
hopping. The lattice is linked to reservoirs through the leads, such a kind of lattice may be 
one dimensional quantum dot array or chains of atoms. In this work we adopt the Keldysh 
NGF formalism rather than other method such as classic or quantum rate equation, since 
the NGF prescription allows one to obtain explicitly analytical current and transmission for- 
mulas, which are valid in both the strong and weak coupling limits. We develop a technique 
to calculate the retarded Green functions for every site based on the equations-of-motion 
method. With the help of the technique we developed, it is very convenient and intuitive to 
calculate the retarded and correlation Green functions for every lattice sites. As the trans- 
mission and current formula is applied to the cases of a single and a double site structure, 
some well-known results are rediscovered, with some new phenomena being observed. The 
case of triple sites shows more complicated and interesting resonant transmission structure. 
If the lattice-lead couplings are asymmetric, the resonant structure of the transmission spec- 
tra for the triple site structure has weak dependence on the arrangement of site levels if they 
are not aligned. Moreover, we derive an analytical expression for the condition for perfect 
transmission of a one dimensional lattice. Such a condition depends on whether the number 
of sites is odd or even, suggesting that the parity feature of the linear conductance of a one 
dimensional system does not arise from the electron-electron correlation effects, since we did 
not consider electron-electron interactions. In the case of symmetric lattice-lead couplings, 
our result is consistent with the parity effect of conductance in monatomic wires based on 
the first-principles calculations Q| and experimental findings Q]. 

we notice that Shangguan et al. studied the differential conductance and charge dis- 
tribution in a linear quantum dot array based on the similar formalism. However, there exist 
some crucial discrepancies between our formalisms and results. First, the prescriptions for 
the calculation of the correlation Green functions such as less-than one are different. After 
repeating the use of the Keldysh formula in the noninteraction case with the help of the 
properties G r G a = (G r — G a )/(2iImG a ), we obtain finally the transmission formula (23), 
which is different formally from the result of Shangguan et al.j^J, Eq. (27) in their paper, 
without using the mentioned property. Such a prescription is more intuitive and facilitates 
greatly the analysis of the dependence of transmission probability on the parameters of the 
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system especially at resonance, as can be seen in the TV-site case. Second, we study the linear 
transport properties such as resonant transmission, linear conductance, while Shangguan et 
al. investigate the nonlinear ones such as differential conductance, charge distribution. We 
also extend the result to conductance quantization in monovalent atomic chains and the 
localization-delocalization transition in a generic one dimensional system, a fundamental 
issue in the condensed-matter physics. 

The rest of the paper is organized as follows. In Sec. II we formulate the transport in 
a one dimensional lattice and derive the formulas for calculating the current, transmission 
and conductance based on the Keldysh NGF method. In Sec. Ill we use the formulas to 
calculate the transmission and conductance spectra for a single cite, a double site and a 
triple site structures, with a detailed analysis of our results. Also we derive an analytical 
expression for the condition for perfect transmission of A-site lattices as the site levels are 
aligned. Concluding remarks are given in Sec. IV. 

II. MODEL AND FORMULATION 

We are interested in the electronic transport properties of a one dimensional lattice with 
N sites. The lattice is connected to the left(right) lead L(R), which can be described by 
a tunneling matrix element. The leads are considered electron reservoirs with a continuum 
of states filled up to their respective Fermi levels /il and /ir at zero temperature. Such a 
device can be described by a single-band tight-binging model. The on-site energy of site i 
is labelled by q. An inter-site coupling with matrix element U accounts for the electron's 
hopping between nearest-neighboring sites i and i + The coupling between the left (right) 
lead L(R) and the first 1 (last N) site is described by a parameter v k L {vkn)- Then a tight- 
binding Hamiltonian for the lattice in the absence of electron-electron interaction takes the 
following form 

N N-l 

H = e k a at a a ka + ddfdi + ^ (Udf +1 di + t*dfd i+1 ) 

k,a=L,R i=l i=l 

+ E(^4l^i + vl L dta kL ) + Y.( v kR a kndN + v* kR d%a kR ) (1) 
k k 

where a k L(a^ L ) ,a k R(al R ) an d di(df) (i — 1,2, ...N) are the annihilation (creation) operators 

for an electron in the left lead, right leads and at site i. In Eq. (1), the first term is the 

Hamiltonian for the leads, the second term is the Hamiltonian for uncoupled N sites, while 
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the third and last two terms describe the coupling between neatest-neighboring sites and 
the couplings between the lattice and the leads. 



Following Ref. 12j, the current through the lattice is the time evolution of the electron 



number Ni(t) = Y^k a kL(t) a kL{t) in the left lead 

i p Op 

J L = -e < N L >= - < {N L ,H} + >= -Re{J2v kL G< k (t,t)}, (2) 

k 

where the less-than Green function Gf k (t, t') = i < a kL (t')di(t) >. With the help of Dyson's 
equation, the less-than Green function Gf k (t,t') can be written as 

G< k (t,t') = J dhvt^G^ttJg^t') + G<(t,t 1 )g a kL (t 1 ,t% (3) 

In Eq. (3), tf) = i < dt(t')di(t) >, G r n {t,t') = -i6(t - t') < {d 1 (t),df(t')} + > are 

the less-than and retarded Green functions of electrons at site 1, while gl^it, t') is the exact 
advanced (less-than) Green function of electrons in the left lead decoupled from the lattice 

gH(t,t') = Tie{±t T t')e-^^ t - t '\ 
9k L (t,t') = ihieu, - p L )e-<**-'*X t -*\ 

Here f a {x) = [exp(x/kBT) + (a = L, R) is the Fermi-Dirac distribution function of the 
Lead. Substituting Eq. (3) into Eq. (2), changing the sum over k into an integral / depi(e) 
(where pz(e) the density of states in the left lead) and introducing a linewidth function 
T^(e) = 2tt I VkL | 2 Pl( 6 ); we get a compact form of current formula in Fourier space |l2j 

Jl = jS ^ T ^){h(t-^)l G 'n(t) -G?i(c)] +G?i(e)} (4) 

Now we turn to calculate the retarded (advanced) Green function G T {i(e) and the less- 
than Green function Gf^e) of site 1. The essential idea to solve the retarded Green function 
is that we consider first site 2 neighboring to site 1 as a part of the right lead, and then 
consider site i + 1 as a part of the right lead when calculating the Green functions of site i. 
This procedure may be regarded as an extension of the well-known recursion or decimation 
method for calculating Green function. |24( From Dyson's equation G = (g^ 1 — S) _1 , G^(e) 
can be written as 

G 1 n (e) = [g[ 1 (e)- 1 -^n(e)}-\ (5) 
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SI(e) = Z k \v kL \ 2 9l L (e)=AL(e)--T L (e), 



where <?n(e) is the retarded Green function of site 1 decoupled with the system, and 

= ^(e) + S^(e), (6) 

1 (7) 
I 2 Gl R (e). (8) 

G2i?(e) in Eq. (8) is the retarded Green function of site 2 decoupled from site 1 but still 
coupled to site 3. It is given by 



G r m (e) 



J 2R 



with S2J?( e ) = l ^2 I G r m (e). Similarly, the retarded Green function G r iR of site % is 



(9) 



(10) 



*i 1 2 ^i+lii( e )- 



Note that for the last site N, its retarded Green function reads as 



GNR( e ) 



1 ~ ^(6) 



(11) 



where £#(e) = S fe | t> fefi | <7fcfl( e ) = ^n( e ) ~ fr,f?( e )- Substituting the Green function G\ R 
into the exprsession of the Green function G\_ lR (i = 2,3, ■ ■ ■ N) recursively, we finally 
obtain an analytical expression for G r 11 {e) 

1 



G r n (e) 



9n(e) 



-i 



A L + - 2 r L 



(12) 



tN- 



9nn{£) 



Once the retarded and advanced Green functions are known (G a = (G r )*), the less- 
than Green function can be evaluated with the help of the Keldysh formula for the present 
problem 



G< = G r S < G a 



Then 



Gfi 



G r -G a 
1/G a - 1/G r 



^11 ^11 ^< 



1/G?! - 1/G[i 



s~i< ^iR ^iR y< 

~ 1/G? R - «' 

r ,< _ U NR U NR v < 

6 



(13) 

(14) 
(15) 
(16) 



where the self energies are given by 

E< = £< + E< R = S fc | v kL | 2 g< L (e)+ \ h | 2 G^(e), (17) 
= i/L(e-e^)ri+|«i| 2 C?^(e), (18) 

= I U I 2 GfeiuW. (19) 
= S fc | w fcL | 2 ^ L (e) = if R (e - n R )T R . (20) 

Since 1/G a - 1/G r = E r - £ a , combining Eqs. (14) - (16) yields 

G < (e) = JA-riTL + Mt-^r, _ 

J- L + 1 2 

where T 2 = i | ti | 2 [G^e) — ^^(e)] = — 2\t\\ 2 ImG r 2R . Substituting the above expression 
into the current formula (4), one obtains the following Landauer-Biittiker-type formula (l^ 

Jl = -f / ^[h(e - n L ) - f R (e - ^^L-IrnCT^e). (22) 

Equation (22) is the central result of this work. This formula for steady transport is valid 
both in the strong and weak coupling limits. It is also applicable to the nonequilibrium 
situation. In Eq. (22), the term 

2F F 

T(e) = -—^rlmG^e), (23) 

is the transmission coefficient for electron tunneling through the one-dimensional lattice. 
Furthermore, the differential conductance Q in the linear regime jMi — fi R — > + can be 
readily derived as 

It is noted that the above procedure to calculate the related Green functions can be 
extended to two and three dimension lattices with each site having multiple levels in the 
presence of electron-electron interactions treated in the mean-field approximation, which 
will be published elsewhere. In addition, such a procedure is very suitable to numerical 
calculation of the spectral density and tunneling quantities in coupled many site systems. 



III. TRANSMISSION AND CONDUCTANCE 



In this section we will investigate coherent transmission and conductance in a one dimen- 
sional lattice, which may be related to a one dimensional quantum dot/well array or a chain 
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of single atoms. Our interest is mainly to find out how parameters determine the number 
and position of resonant transmission and conductance peaks. 



A. Single site 

In the case of single site, the Green function and self-energy take the following simple 
form 

G r n (e) = (e-e 1 -A+~T)-\ 

r 2 = T R , (25) 

where A = + Ar, T = Tl + Tr. Then the transmission and conductance are [l^ 



T(e) = r L r^( e 



e 2 r L r Rr d 



m —Ref^ + A-efiL + iT/2) 

^^'i^^ 1 (26) 



where ^i(e) = [(e -e x - A) 2 + (T L + Tr^/A]' 1 , and is the trigamma function [25(. It 
is obvious that the transmission coefficient is of the Breit-Wigner type. Only if the coupling 
between the site and the two leads is symmetric, i.e., T L = T R , perfect resonant transmission 
(T = 1) can occur at the renormalized dotsite level e\ + A, while for asymmetric coupling 
Tl 7^ Tfi; the transmission coefficient T is always less than 1. The larger the asymmetry, 
the smaller the transmission coefficient. Thus conductance quantization can be achieved for 
the 1-site system in the symmetric lattice-lead coupling case. 



B. Double sites 



In the case of two coupled sites, it is expected that the competition between the coupling 
of the sites and the sites and leads play a crucial role on the transmission and conductance. 
In the absence of coupling to the leads, the level associated with two equal sites will split 
into two levels due to the coupling between them. The separation between these two split 
levels is proportional to the coupling strength. When the sites are connected to the leads, 
one should consider the sites and the leads as a single system, with the result that now the 
size of the level splitting depends on both the couplings between sites and leads. 
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Consider first a double site structure without coupling to the leads. Assuming the cou- 
pling strength between these two sites is t, the hamiltonian becomes €\d\ d\ + e 2 (i 2 <i 2 + 
tdid,2 + i*c?2~(ii. Then the retarded Green functions for the two sites are 

1 



G\{e) 



e-ei-iO+- |t| 2 /(e-e 2 -iO+)' 
1 



" e-e 2 -iO+-\t\*/(e-e 1 -iO+Y (27) 
One can readily find that the above two retarded Green functions have the same two poles 
at 

e 1 + e 2 ± v /(ei-e 2 ) 2 + 4|t| 2 
e = > |_L_ (2g) 

Note that when e\ = e 2 , the level separation is simply 2\t\. The above analysis clearly shows 
that the larger the coupling, the bigger the separation between the two split levels. When 
the sites are connected to the two leads, this new coupling may modify the effective coupling 
between the sites as will be discussed in what follows. 

The associated retarded Green functions in the presence of coupling with the two leads 
read 

G r n (e) = [e-e 1 -A L + -Y L - \t\ 2 G r 2R (e)}-\ (29) 
G r 2R (e) = (e-e2-A R + l -T R )-\ (30) 

and 

T 2 = r*|*| 2 /B(e), (31) 

where BU) = (e - e 2 - A R ) 2 + T 2 R /A. Substituting Eqs. (29)-(31) into the expression (23) 
yields J26J 

T(e) = r L r R \t\ 2 F 2 (e) (32) 

where 

13(c) 

^ = [(e - ei - A L )B(e) - \t\\e - e 2 - A R )] 2 + \T L B(e) + |t| 2 r fi ] 2 /4' (33) 
We now consider the case when the levels of the isolated sites 1 and 2 are the same 
(ei = e 2 = eo). In the case of symmetric coupling, i.e., Tl = T R = T, Al = A R = A, the 
condition for perfect resonant transmission is that 

(e-eo- A) 2 + T 2 /4= \t\ 2 (34) 



has real roots. This equation also determines the number and position of the resonant 
transmission peaks. Obviously, the condition for perfect resonant transmission is \t\ 2 > 
T 2 /A. There will be then only one perfect resonant transmission peak located at e = e + A 
when \t\ 2 = T 2 /4. If \t\ 2 < T 2 /4 , there is just one imperfect vresonant peak (T < 1) 
pinned at e = eo + A. In the case of \t\ 2 > T 2 /4, two perfect transmission peaks exist at 
e± = e + A ± ^\t\ 2 - T 2 /4. These features can be clearly seen in Fig. 1 (a), where we can 
also appreciate that the lineshape of all the transmission peaks is Lorenzian. Here and in 
all figures following energies are in arbitrary units. 

Next we consider what would happen if the coupling between the sites and the two leads 
becomes asymmetric. In this case, one can find from the expression of the transmission 
coefficient (32) that the condition for perfect transmission is \t\ 2 = TlTji/4 and only one 
perfect transmission peak can be expected at e = e + A. In all other cases, the transmission 
coefficient is less than I. In Fig. 1 (b) we show the transmission coefficient for various inter- 
site coupling constants t for T L = 0.25 and T R = 4. If the level shift A L and A R induced 
by the left and right lead are not the same, the lineshape of the transmission peaks is non 
Lorenzian, as is apparent in the case A^ = 0.025, Ar = 0.4 also shown in the figure. Figure 
1 (c) and (d) are the plots of the transmission coefficient for increasing inter-site coupling t 
in the symmetric and asymmetric coupling cases, respectively. When the levels of the sites 
are not aligned. It is then expected that no perfect resonant transmission exists. One can 
observe that two symmetric transmission peaks always resolved in the symmetric coupling 
case. As the inter-site coupling t increases, the maximum value as well as the separation 
of the transmission peaks increases. For asymmetric coupling, one transmission peak can 
be expected when the coupling between the sites is weak, and two asymmetric transmission 
peaks when it is strong. In addition, with increasing interdot coupling t, the value of the 
transmission coefficient increases first and then decreases after it is saturated. The position 
of the principal transmission peak is closer to the level of the site with smaller coupling to 
the lead. 

This complex behavior reflects the competition between two resonances derived from the 
ground state of each site when far apart. As they come together and when coupling increases, 
the resonances approach the ground and first excited states of the compound system, which 
in general separated in energy. If the coupling is weak, on the other hand, the peaks will 
only be resolved if their width is smaller than the energy difference Ae = |e 2 — ei|. If the 
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couplings to the leads is asymmetric the competing resonances will be dominated by that 
whose associated wavefunction is concentrated on the site region must weakly coupled to 
the lead. 

To sum up, we have investigated in detail the resonant structure of the transmission 
coefficient of a coupled double site structure in the linear regime. We found that when 
the levels of the two sites are aligned, the condition for perfect resonant transmission is 
\t\ 2 > T 2 /4 in the case of symmetric site-lead coupling, and |t| 2 = r^r #/4 in the asymmetric 
case, which is consistent with the derivations of Ref. 2^|. Once there exists a mismatch 
between the two levels of the sites, no perfect resonant transmission peak can be expected. 
In the case of symmetric site-lead coupling, the value of the transmission coefficient increases 
with increasing inter-site coupling, until it is saturated. For asymmetric site-lead coupling, 
the transmission coefficient increases first and then decreases after it reaches a maximum, 
as the inter-site coupling is increased. Moreover, in the asymmetric case, the splitting of the 
transmission peak requires a stronger inter-site coupling t. Asymmetry between the two split 
transmission peaks can be observed when the two levels are not aligned for asymmetric site- 
lead coupling, provided the inter-site coupling is strong enough. In addition, different level 
shifts of the two sites induced by the coupling to the leads will introduce a non Lorenzian 
lineshape of the transmission peak(s). The results imply that no perfect transmission and 
then no conductance quantization can be expected for double site systems in the symmetric 
coupling case (inversion symmetry). 



C. Triple Sites 

A coupled triple site structure is expected to contain richer physics than a coupled double 
site system, since it permits more interesting arrangement of the energy levels in each site and 
the competition between inter-site couplings. Surprsingly, electronic transport through the 
triple sites has been less investigated in the past both theoretically and experimentally [28j]. In 
this subsection, we study the detailed dependence of the number and profile of transmission 
peaks, on the parameters of such a structure. 

For an isolated triple site system, the hamiltonian can be written as 

H = tidfdi + t 2 d 2 d 2 + e 3 djd 3 + t L d\d 2 + t* L d\d x + t R d^d 3 + t* R d^d 2 , (35) 
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where ti and t R are the couplings between sites 1 and 2, and 2 and 3, respectively. The 
related retarded Green functions are readily derived as 

1 



G[(e) = 



e — ei — i0+ 



I*lI 2 



-e 3 -i0+ 

G5(c) = ! 



e - e 2 - 20+ - |i L |2/(e - Cl - *0+) - \t R \ 2 /(e - e 3 - *0+) ' 
1 

e — e 3 — i0+ — 



G5(c) = M (36) 



e— 62— i0+- 



-i0+ 



when the levels of three sites are aligned, i.e., e\ = e 2 = e 3 = eo, each of these three 
retarded Green functions has the same three poles at e = e , e = e + \J\th\ 2 + |£r| 2 and 
e = e — ^I^lI 2 + I^rI 2 - When the three electron levels are mismatched, there exist in 
general three poles in the above retarded Green functions, that is to say, the position of the 
original levels are moved due to the existence of the inter-site couplings. As the couplings 
between the structure and the leads are turned on, the Green functions for calculating the 
transmission coefficient become 

GIi(c) = [e - d - A L + -Y L - \t L \ 2 G r 2R (e)}-\ (37) 

G r 2R {e) = {e-e 2 -\t R \ 2 Gl R )-\ (38) 

G r 3R (e) = (e-e 3 -A R + l -r R )-\ (39) 

and 

T 2 = -2\t L \ 2 ImG r 2R . (40) 
Then the transmission coefficient becomes 

T(e) = T L T R \t L \ 2 \t R \ 2 ^ 3 (e), (41) 

where 

T (A = ^ (AO) 

^ > C (e)[(e -e 1 -A L - \t L \ 2 ReG r m y + (T L - 2\t L \^ImG r 2R f /4] 1 ] 



with 



B(e) = (c - c s - A*) 2 + ry4, 

C(e) = [(e-e 2 )B(e)-\t R \ 2 (e-e 3 -A R )] 2 +\t R \'T 2 R /A, 
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and 

ReG r 2R (e) = [(e-e 2 )B(e)-\t R \ 2 (e-e 3 -A R )] B{<) 



C(e)' 



ImG 2R (e) = — ^-y. (43) 

Now we consider the case when the three levels of sites are aligned and neglect the energy 
shifts Al and A R . At resonance e = eo, one has 

B(e) = I* /4, 
C(e) = \t R \*T* R /4, 
ImGr 2R (e) = (44) 

then one can easily derive the condition for perfect transmission 

£ i 2 = («) 

tR 1 R 

Further analytical results in this case are cumbersome, so we shall next provide some numer- 
ical results on the transmission spectra of the coupled triple site structure in various cases. 
Figures 2 and 3 show the transmission spectra for symmetric coupling (r^ = T R = 1), 
and asymmetric coupling (T L = l,T R = 4) to the left and right leads, respectively. In 
the figures four kinds of arrangements of energy levels are presented: (a) aligned levels 
ei = ^2 = £3 = 10, (b) ladder levels e\ = 9, e 2 = 10, e 3 = 11, (c) y-type levels e\ = 9, 
e 2 = 11, e 3 = 9 and (d) two aligned neighboring levels e x = 9, e 2 = 9, £3 = 11. From 
Fig. 2 (a), one finds that there exists only one perfect transmission peak located at e = eo 
when the inter-site couplings are equal and small. As ti = t R are increased, the perfect 
transmission peak splits symmetrically. Once the energy shifts are taken into consideration, 
the transmission peaks are no longer equally spaced (dotted line in Fig. 2 (a)). Since Eq. 
(45) is no longer satisfied, there is no perfect transmission when 7^ t R . 

In the case of levels arranged in a ladder sequence (Fig. 2 (b)), one perfect transmission 
peak with two low shoulders is obtained for ti, — t R — 0.5. As they both reach the value 2, 
the shoulders become two imperfect transmission peaks and the perfect transmission peak 
at the center is widened. When the inter-site couplings are not equal, the transmission peaks 
are suppressed and arranged into a ladder type. When the three site levels are in the l^-type 
disposition, interesting transmission spectra are observed. If the inter-site couplings are 
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equal and small, one sharp perfect transmission peak and one broad imperfect transmission 
peak are seen. As the inter-site couplings increase, the imperfect transmission peak splits 
and the resonances move further apart. Similarly, the inequality of t L and t R decreases the 
transmission in all three peaks. In fact the value of the transmission coefficient through split 
electron levels is determined by the extension to the leads of the wavef unctions for these split 
electron levels. Then the above phenomena can be similarly explained as we did in the case 
of two coupled sites. Figure 2 (d) shows the transmission spectra for e\ = 9, e 2 = 9, e 3 = 11. 
One can see two asymmetric imperfect transmission peaks as = t R = 0.5. When they 
reach the value 2, three imperfect transmission peaks of different height are discerned. As 
in the case of three ladder levels, the asymmetry between two inter-site couplings and tR 
rearranges these three transmission peaks into a ladder. Notice that no perfect transmission 
peak even appears in this case. 

Now we investigate how the transmission spectra are modified when the site-lead cou- 
plings are asymmetric. One can find from Fig. 3 that perfect transmission exists only in 
the case of three aligned levels under the condition for perfect transmission. Figure 3 (a) is 
the transmission spectra for different inter-site couplings when the levels of three sites are 
aligned. One imperfect transmission peak is seen as tL — tR and one perfect transmission 
peak along with two narrow shoulders is observed as tL — 1, tR — 2, which satisfies the con- 
dition for perfect transmission (45). Also one can find different resonant structures for the 
cases \t L /t R \ 2 > T L /T R and |£l/£r| 2 < T L /T R . When the three levels are not aligned, no per- 
fect transmission occurs in the case of asymmetric site-lead couplings, which is demonstrated 
in Fig. 3 (b)-(d). The transmission spectra in the case of levels arranged in a ladder exhibit 
one imperfect transmission peak and one shoulder structure as = t R = 0.5. Increasing the 
inter-site couplings, the shoulder becomes gradually distinguished as a transmission peak. 
The competition between t L and t R strongly modifies the resonant structure of the trans- 
mission spectra. Comparing Fig. 3 (b)-(d), one may find that the resonant structure of the 
transmission spectra are similar and mainly display two transmission peaks structure. It 
suggests a fact that, as the site- lead couplings are not symmetric, the resonant structure of 
the transmission spectra has weak dependence on the arrangement of the electron levels of 
sites if they are not aligned. 

From the above analysis, the transmission spectra of three coupled sites is more compli- 
cated and contains more physics than that of two coupled sites. The main features of the 
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transmission spectra for a triple site structure are as follows: (1) the transmission spectra 
may contain just one, two or three resonant peaks, (2) the resonant structure of the trans- 
mission spectra is strongly dependent on the arrangement of the electronic levels when the 
site-lead couplings are symmetric, while weak dependence on the arrangement of the elec- 
tron levels if they are not, (3) the resonant structure of the transmission spectrum strongly 
depends on the symmetry between inter-site couplings, (4) perfect transmission and thus 
conductance quantization can be achieved for a triple site system with inversion symmetry 
when coupled symmetrically to reservoirs. These conclusions hold for the linear conductance 
Q as well, since the conductance and the transmission coefficient possess the same resonant 
structure at low temperatures. Our theoretical results are consistent with the experiment 



observations in the conductance of triple sites 



D. N Sites 



28). 



As can be found from our previous analysis, the condition for perfect transmission is 
different for one, two and three site structures. It gives us a hint that the condition for 
perfect transmission would depend on whether the number of sites in a ID lattice is even 
or odd. In what follows, we derive the condition under which the perfect transmission 
happens when the electron levels of the lattice are aligned and energy shifts are ignored, 
i.e., e\ — 62 — ■■■ — ejsr — eo, and = A R = 0. From the expression for the transmission 
coefficient (23) and for T 2 , one has 



4r L 
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to 


ImG r 2R 


r L - 


2| 


k 


2 ImG r 2R 



(46) 



-2r L 
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2 ImG r 2R 


(e - e - 
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2 ReG r 2R ) 2 + (T L -2\t l 


2 ImG 2R ) 2 / A 



At resonance e = €q, the real part of all the retarded Green functions becomes zero. Then 
the transmission coefficient is given by 



(T L -2\t 1 \ 2 ImG r 2R ) 2 1 ) 

Perfect transmission through the lattice, T = 1, is obtained if 

T L = -2\ti\ 2 ImG r 2R . (48) 
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Notice that since g\ { = {% = 2, 3, • • • TV) at resonance e = eo, one can obtain the following 
expression for ImG r 2R 

ImG r R = _^R ^---t N -2 l2 . {4Q) 
2 t 2 t^ ■ ■ ■ tw— 1 

JmGL, = -A| Mi^zl | 2; ive,en. (50) 
Then the condition for perfect transmission is 

N odd (51) 

N even. (52) 
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If the interdot couplings are the same, the condition becomes 

T L = T R ; N odd (53) 

\ tN \ 2 = TV even. (54) 

t 4 

Thus one can conclude that for a one dimensional lattice, the condition for perfect trans- 
mission is dependent on the parity of the sites of the lattice, i.e., whether the number of sites 
is odd or even. Equation (51) suggests that for a chain with an odd number of sites and in- 
version symmetry (i.e., T L = T R = T,t\ = t/v-i,t2 = tjv-2; etc.) perfect transmission will be 
automatically satisfied at the middle of the band or level group. Thus the linear conductance 
is quantized to the value 2e 2 /h. This is not the case when N is even. From the Eqs. (52) 
and (46), we find a transmission coefficient An/ (1 +n) 2 with n = I rt 2 tT'fiv~2 ^ ^ ess ^ nan um ty 
and a conductance smaller than 2e 2 /h. The even-odd feature appears in transmission and 
conductance in the absence of any electron interactions, thus the parity feature is not due to 
many particle effects. Our argument also proves that when the system is symmetric under 
inversion, the state at the middle of the band is always delocalized, regardless of the amount 
of disorder that respects such symmetry condition. Eqs. (51) and (52) defines a broad class 
of correlations in the disorder yielding derealization of ID disordered systems. This parity 
effect of conductance is consistent with the predictions in the monatomic wires based on the 
first-principles calculation^] and confirmed by the experimental observations [l^ . 

The dependence of transmission or conductance on the number of sites can be explained 
qualitatively as follows. To simplify our discussion, we constrain ourselves to the case when 
the inter-site couplings are the same, which corresponds to the monovalent atomic wire 
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easeful]. When the number N of sites is odd, one finds that the current through the middle 
site (N + l)/2 can be written as 

Jo = / TpLMe - n' L ) - f R {e - n'n)} ^ 10 ^* ImGVi^fe), (55) 

h J 27T 1 LD + 1 RD 2 2 

where T LD = —2t 2 ImG r N ^ 1 T , T RD = —2t 2 ImG r N+3 D , with G n-i t and Gn+3 r are the Green's 

— 2 — — 2 — 2 2 

functions decoupled from the middle site (N + l)/2. In Eq. (55) we denote the chemical 
potential of the site left (right) to the middle site by /i' L (fjf R ) which equals to the chemical 
potential hl (i^r) in the equilibrium. The conductance at resonance can then be expressed 
as 

e ° = -(T LD+ r RD r (56 » 

Obviously, the conductance reaches ^ only if T LD = T RD , which implies T L = T R from the 
expressions of Tld and T R d. When the site number N is even, one can divide the lattice 
into two parts and then the current from the site N/2 to (N + 2)/2 can be written as 2^ 

o P r A f At 2 ImG r N r ImG r N+2 R 



* L ^E±l R \ 



The conductance at resonance (ReG r = 0) is 

2e 2 4:t 2 ImG r N L ImG r N+ 2 R 



h 1 1 + t 2 ImG r N T ImG r N+2 D 



2 ' 



(58) 



If Qe = one needs t 2 = (ImG r N ImG r N+2 n ) 1 , which eventually yields t 2 = 

ft — L ~^2~ H 

after simple calculation. It is in fact the perfect transmission condition in the case of double 
sites. Physically the even-odd dependence of conductance reflects the different conditions 
of constructive interference in the different circumstances. For the lattice with odd number 
of sites, one can view the array as a single site coupled with two renormalized leads. If the 
number of sites is even, the lattice can be considered as a two renormalized normal metal 
contact or a double sites coupled with two renormalized leads. The condition for perfect 
transmission is more strict in the later case. 

Finally, we would like discuss possible effects arising from electron-electron interaction. 
It is expected that the main results may hold as well. In the presence of on-site Coulomb 
interaction, an additional term of self-energy will be introduced for the Green function of each 
site. The effect is just shifting the site level and splitting the resonance position at absolute 
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zero temperature, since elastic scattering can only be expected at zero temperature. [3(J The 
central formula (23) remains formally the same in the presence of interactions as long as the 
ground-state of the system possess Fermi liquid properties. 



IV. CONCLUSIONS 



In summary, using the Keldysh nonequilibrium Green function method, we have derived 
the formulas to calculate the transmission coefficient, current and conductance of a chian of 
coupled N site system. An effective and convenient procedure to calculate recursively the 
retarded(advanced) and lesser (greater) Green functions has been developed. Based on the 
formulation developed, we have analyzed the transmission spectra of just single site, double 
site and triple site structures in detail, obtaining some well-know results and finding some 
new features in the transmission spectra of double and triple site systems. When the electron 
levels of N sites are aligned, we have obtained an analytical expression for the condition 
for perfect transmission, demonstrating the even-odd parity effect in the transmission of a 
generic one dimensional lattices with inversion symmetry. 
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FIG. 3: Transmission spectra of a triple site structure in the case of asymmetric site-lead couplings 
(r L = l,T R = 4): (a)ei = e 2 = e 3 = 10; (b)ei = 9,e 2 = 10, e 3 = 11; (c)ei = 9,e 2 = 11, e 3 = 9 and 
(d)ei = 9,e 2 = 9,e 3 = H- 
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fig. 3 Zeng et al. 
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